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Abstract 

We study the complete one-loop electroweak(EW) corrections to the processes of 
single charged Higgs boson production associated with a neutral Higgs boson(/i°, H°, A°) 
and a gauge boson in the framework of the minimal supersymmetric standard 
model(MSSM). Numerical results at the SPSla' benchmark point as proposed in 
the SPA project, are presented for demonstration. Wc find that for the process 
e+e~ — > h'^H^W^ the EW relative correction can be either positive or negative and 
in the range of —15% ~ 20% in our chosen parameter space. While for the processes 
e+e~ — > H'^{A^)H^W^ the corrections generally reduce the Born cross sections and 
the EW relative corrections are typically of order —10% ~ —20%. 
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1 Intorduction 



Even though the Standard Model(SM)[Tl [2] has been successful in describing the strong, 
weak and electromagnetic interaction phenomena at the energy scale up to 10^ GeV, it 
has been known for a long time that it is still considered as a low-energy effective theory 
suffering from a number of theoretical difficulties, such as hierarchy problem. It is likely 
that at the higher energy scale the minimal supersymmetric standard model (MSSM) [31 [U 
[S] is the most attractive candidate among various extensions of the SM, which is expected 
to resolve theoretical difficulties existing in the SM. By adopting two Higgs doublets to 
preserve the supersymmetry in the MSSM, five Higgs bosons are predicted. They are two 
CP-even Higgs scalar bosons and H^, one CP-odd Higgs boson and a charged Higgs 
pair 

The charged Higgs boson plays an important role in the extensions of the SM, its 
discovery in experiment will unequivocally imply the existence of the physics beyond 
the SM. If the charged Higgs boson really exists in nature then it is mostly possible 
to be discovered in the coming few years at the CERN Large Hadron collider (LHC) due 
to its TeV scale colliding energy, but the further accurate measurements of masses of 
Higgs bosons and investigations of the properties of Higgs bosons [6J, which is crucial for 
verification of electroweak (EW) symmetry breaking mechanism, need be best performed in 
a clean environment of e^e~ collisions, such as at the future International Linear Collider 
(ILC) [7]. The study would be of the utmost importance for establishment of the MSSM 
or some of more general supersymmetric models as the theory of the EW interactions. In 
addition, further study of the complete Higgs spectrum is the cornerstone of the definitive 
test of mechanism of electroweak symmetry breaking. In discovery and detailed study of 
the properties for these Higgs bosons, there are some important processes for neutral and 
charged Higgs bosons at linear colliders, such as, e^e^ Z* ^ h?{H^)Z^ , e~^e^ ^ Z* — > 
h°{H°)A^, e+e" ^ iyi?W+*W~* iyi?h°{H^), e+e" ^ H+R- and so on[8j. The lightest 
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neutral Higgs boson in the MSSM has the mass upper hmit of about 140 GeV, and it 
decays dominantly into bb pairs. For the heavier neutral Higgs bosons as a result 

of the strong enhancement of the coupling between Higgs boson and down-type fermions, 
both and will decay almost exclusively into bb (~ 90%) and r^r^ (~ 10%) pairs 
(The decays H'^,A'^ — > tt are possible only when they are kinematically allowed), and on 
the contrary all other decays are strongly suppressed^ [ini [HI 1121 |T3]. For the mass of 
charged Higgs boson being blow 170 GeV — ISOGeV (i.e., the H~ — > tb decay threshold), 
the dominant decay mode of the charged Higgs boson is t^Vt- But when the mass 

of charged Higgs boson is above the tb decay threshold, the charged boson 

mainly decays into tb pair leading to Wbb final states 

If H^H^ pair production at e~^e~ linear colliders(LC) is kinematically allowed, the 
dominant production process of charged Higgs bosons is the pair production via exchang- 
ing a virtual photon and Z boson e~^e~ — > 7*, Z* — > H~^H~[8\. For a fixed e^e~ colliding 
energy, its cross section is related only with the charged Higgs boson mass. Another im- 
portant production mode is the production of single charged Higgs boson, especially when 
the charged Higgs boson is too heavy to be produced in pair at linear colliders. Among 
various single charged Higgs boson production channels, the bottom quark associated 
production process of e~^e~ — > bbH^W^[14\ is experimentally interesting since it can be 
yielded by many topologies, such as the tt and H^H^ pair productions followed by their 
corresponding sequential decays. The e^e~ bbH^W^ process can be also induced from 
following three single charged Higgs boson production processes ([1]) and the sequential 
decays h^{H^,A^) bb, 

e+e' ^ ^^H^W^, ($° = H^, A^). (1) 

S. Kanemura, S. Moretti, et al, studied the processes of single charged Higgs boson 
production channels in the MSSM at LCs, including the processes ([TJ at tree-level in 
Ref . [14t I15j . It was found that there are regions of parameter space where it is acces- 
sible beyond the kinematic limit for pair production, and these processes are excellent 

3 



channels in studying the interactions among neutral Higgs boson(/i'^, or ^4"), charged 
Higgs boson and gauge boson(M/^^), and is a useful alternative channel for studying the 
phenomenology of charged Higgs boson. 

In present work, we will study the EW corrections to the processes ([1]) including 
complete one-loop diagrams in the MSSM. The paper is organized as follows: In Sect. 2, 
we present the analytical calculation of the complete one-loop EW radiative corrections to 
e+e" ^'^H^W'^i^'^ = hP,H^,A^) processes in the MSSM as well as the treatment of 
charged Higgs boson resonance. The numerical results and discussions are given in Sect. 3. 
Finally, we give a short summary. 

2 Calculations 

2.1 Conventions and notations in analytical calculations 

We denote the processes as 

e+(m) + e-fe) - ^°(P3) + ^+^4) + W-{p,), (cDO = /lO, ^O). (2) 
The lowest order cross section reads 

a,.ee = J^/E|-^*-|'^^3 (3) 

' ^ spin 

where pi is the c.m.s. kinematical momentum of the incoming positron. The integration 
is performed over the three-body phase space of final particles <^^H~^W~ , the bar over 
summation indicates averaging over initial spins. The phase-space element is defined 

by 

*=*«(«+p.-|:p,)np|i^ (4) 

In our calculation, we use dimensional reduction(DR) regularization scheme, which is 
supersymmetric invariant at least one-loop level, to regularize UV divergences and adopt 
on-shell(OS) conditions (neglecting the finite widths of the particle) to renormalize the 



relevant fields [16]. We adopt the 't Hooft-Feynman gauge and the definitions of one- loop 
integral functions in Ref.[17j in the calculation of one-loop diagrams. The Feynman dia- 
grams and their corresponding amplitudes are created by FeynArts 3 [18j automatically, 
and the Feynman amplitudes are subsequently reduced by F0RM[19\. We depict the 
lowest-order Feynman diagrams in Fig.l. The numerical calculation of the two-, three- 
and four-point integral functions are done by using FF package pOj. The implementations 
of the scalar and the tensor five-point integrals are carried out exactly by using the Fortran 
programs as used in our previous works [2T] with the approach presented in Ref. |22| . Due 
to the fact that the Yukawa coupling of Higgs / Goldstone to the fermions is proportional 
to the mass of the fermion, we ignore the contributions of the Feynman diagrams which 
involve the Yukawa coupling between any Higgs/ Goldstone boson and electrons. 




Figure 1: The tree-level Feynman diagrams of the processes e~^e ^^H~^W = 
The 0{aew) virtual correction to the cross section can be expressed as 

(^virtual = (^treeS virtual = „, ^ , ^ / '2^*^3^-Re(A^(ree-A^|^jrt«a«) (^) 

spin 

where virtual is the renormalized amplitude involving all the one- loop electroweak Feyn- 
man diagrams and corresponding counter-terms p3] . 
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The definitions and the exphcit expressions of these renormalization constants can be 
found in Refs. [Mll25t l26|. Here we just hst the definitions, which do not appear in Ref. 
as below: 
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■ Re 
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5^0 ^0 



(6) 



The Higgs tadpole parameters fty^o^o, hQo^o, ft^/o/^o and bjjOHO are defined and ex- 

pressed as in Ref. [25]. The operator Re takes only the real part of the loop integrals and 
does not affect the possible complex couplings. 

2.2 Disposal of the charged Higgs boson resonance 




(a) 



(b) 



Figure 2: Feynman diagrams with charged Higgs boson resonance for process e~^e 
h°H+W- 

The disposal of resonance in the calculation of the EW correction is similar with 
the method described in Ref.|27j. For the process e'^e~ — > h^H^W~ in TeV scale colliding 
energy, the amplitude of the Feynman diagram with internal line of (shown in Fig. 2(a), 



denoted as ^^ireeO^H- = 0)) suffers from resonance effect when S34 = (pa + P4) 



M' 



H- 



due to the real charged Higgs boson decay H~ — > h^W~. We have to introduce the 
complex mass to do the following replacement in the resonance propagator of the 
amplitude Mj^^^iTfj- = 0), 

1 1 



where Tjj~ is H~ decay width. Then we get the amplitude at the tree-level as 



(7) 



^^r'ee = ''2 r •^^e"e(r//- = 0). (8) 

The modified amplitude M^ree is safe amplitude being free of the H~ resonance sin- 
gularity. We know that at the tree-level calculation this replacement does not induce the 
gauge invariant and double-counting problems[27], but we should remember that after this 
treatment, there includes the contribution from the imaginary part of charged Higgs boson 
self-energy(S^ ^ ) at higher order. In order to avoid the double-counting in calculation 
up to one-loop level, we should subtract the higher order contribution from the tree-level 
amplitude A^t^ee- practical calculation, we separate some special Feynman dia- 

grams(shown in Fig. 2(b)) from entirety. Each of these diagrams in Fig. 2(b) contains a 
one-loop of charged Higgs boson self-energy, or the corresponding counterterm. Then we 
can express the one-loop level amplitude as below: 



M^Z-ZoH" = [Mone-loopiTH- = 0) - A^^-^- (F^- = 0)] _^ 
" ]^"^"(S34) - [S^"^"(M^_) -ft^o^o 



S34 - M2 



+ 6Z 



H~H- 



H- 



M?~e, (9) 



where the first term on right-handed side is just the amplitude at one-loop level af- 
ter subtracting the the contributions of Fig. 2(b). The one-loop(or higher) order con- 
tributions are included by doing the replacement shown in Eq.Q. The real part of 
^ {Mjj_) — bj^oj^o is just equal to 6M'^_, the mass counterterm of charged Higgs 



boson, while Ti^ ^ (.-^h- ) ~ ^a^a" receives its imaginary part in the resonance from the 
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subtraction of the iMjj-T jj- contribution already contained in Mfree- calcula- 
tion, charged Higgs boson width was evaluated using the HDECAY package p8]. and UV 
divergence cancelled numerically. 

2.3 Real photonic corrections 

The 0{aeui) virtual correction suffers from both ultraviolet (UV) and infrared(IR) diver- 
gences. It is checked numerically that the UV divergence vanished after renormalization. 
IR divergence comes from virtual photonic loop correction. It can be cancelled by the 
contribution of the real photon emission processes. We denote the real photon emission 

as 

e+(R)+e-(p2)^^°(p3)+^+M + ^-(p5)+7(M' (^° = /i°,^°,^°). (10) 

We adopt the general phase-space-slicing method [25] to separate the soft photon emis- 
sion singularity from the real photon emission processes. By using this method, the 
bremsstrahlung phase space is divided into singular and non-singular regions by the soft 
photon cutoff /S.E^ , then the relative correction of the real photon emission is broken down 
into corresponding soft term 5soft with radiated photon energy < AE^ and hard terms 
Shard with > AE^, where k^ = \J\k^\'^ + and is the mass of photon, which is used 
to regulate the IR divergences existing in the soft term. The IR divergence from the soft 
contribution cancels that from the virtual correction exact ly [30l 131 j . Therefore, the sum of 
the virtual and soft cross section is independent of the IR regulator m^. In the calculation, 
The phase space integration of the processes e^e~ <^^H'^W~^, {(^^ = h^,H^,A^) with 
hard photon emission is performed by using the program GRACE [32] . The hard photon 
emission cross section ahard is UV and IR finite with ky > AE^y. Finally, we can get the 
UV and IR finite 0{aew) correction. 



8 



3 Numerical Results 



In the numerical calculation, we use the following SM parameters [23] 

me = 0.510998902 MeV, = 105.658369 MeV, rrir = 1776.99 MeV, 

rriu = 66 MeV, rUc = 1.2 GeV, mt = 178 GeV, 

ma = 66 MeV, = 150 MeV, mj, = 4.7 GeV, 

mw = 80.425 GeV, mz = 91.1876 GeV, OewiO) = 1/137.036. (11) 

Here we take the effective values of the light quark masses (rriu and m^) which can repro- 
duce the hadronic contribution to the shift in the fine structure constant Q^e^u(^T^^) [33]. As 
a numerical demonstration, we generally adopt the relevant supersymmetric parameters 
from the CP and R-parity invariant MSSM reference point SPSla' which is defined in the 
SPA project [34l I35j . This point satisfies all the available high- and low-energy precision 
data and both the bounds for the masses of the SUSY particles and the bounds from cos- 
mology. As we use the on-shell renormalization scheme, we have to transform the SPSla' 
DR input parameters V into on-shell parameters P*-*^ . This transformation is simply per- 
formed by subtracting the corresponding counterterms i. e. "P*-*^ = V{Q) — 5V{Q). The 
input parameters for the reference point SPSla' are defined in the DR scheme at the scale 
(5 = 1 TeV , which are a set of mSUGRA parameters listed below: 

Ml/2 = 250 GeV, Mq = 70 GeV, Aq = -300 GeV, (12) 
sign{p) = +1, tan /3(M = 1 TeV) = 10. (13) 

The relevant parameters of SPSla' point are listed in Table [TJ 

Since the difference between DR and OS mass values for sparticle and Higgs bosons is 
of high order, either of them can be used. In our paper we generally adopt the DR masses 
of sparticle and Higgs boson, but when we study the cross sections as the functions of M40 
or take M^o / 424.9 GeV, we still adopt the SPSla' MSSM spectrum and parameters 
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Mass(GeV) 
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400.5 


hi 


506.3 


xl 


413.9 




545.7 


xt 


183.7 






xt 


415.4 







Table 1: Mass spectrum of supersymmetric particles and Higgs bosons in the reference 
point SPSla'. 

except Higgs sector. In the later case the masses of the MSSM Higgs bosons are fixed by 
taking into account the significant radiative corrections in the MSSM by using FormCalc 
package with the input parameters: M^o and OS tan/3 value(i.e., tan/3 = 10.31 at SPSla' 
point) [36]. We checked in our calculation that the total cross section is independent on the 
IR regulator parameter and the soft cutoff 6s- In following we take = 10^^ GeV 
and 5s = 10"^. 

The Born and 0{aeu]) electroweak corrected cross sections at the SPSla' point for 
process e'^e~ — > h^H^W~ as the functions of c.m.s energy in the energy range from 
900 GeV to 1.8 TeV are illustrated in Fig. 3(a). The tree-level total cross section is drawn 
in solid line, while the full one-loop EW corrected cross section in dashed line. We can read 
from the figure that at the position of ^/s = 1.3 TeV, the Born and EW corrected cross 
sections are 0.019 fb and 0.015 fb, respectively. The corresponding relative corrections 
are depicted in Fig. 3(b) and we can see that the absolute relative correction is very large 
when the colliding energy goes close to the vicinity of the threshold for process e~^e~ — > 
h°H+W~. Fig.3(b) shows that when ^/s = 900 GeV, 6totai is about -36%, but when the 
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-y/s goes to 1.3 TeV, the absolute relative correction declines to the value of 19.2%. 




Figure 3: The Born and full 0{aeu)) electroweak corrected cross sections as well as the 
corresponding relative corrections Stotal as the functions of the c.m.s energy -y/s at SPSla' 
point for process e'^e~ — h^H'^W'~. (a) cross section versus s/s. (b) Stotal versus -s/s. 

In Fig.4(a), we present the Born cross section atree and one-loop electroweak corrected 
cross section atotal with M^o = 200 GeV and related MSSM parameters being taken from 
SPSla' scenario except Higgs sector for reaction e+e~ — ^ hPH'^W~. In this case the 
mass of the lightest neutral Higgs boson Mf^o is equal to 117.2 GeV, which is fixed in 
the MSSM by using FormCalc program with the input parameters: M^o = 200 GeV and 
OS tan/3 = 10.31. In contrast with the plots in Fig.3(a-b), the complete EW one-loop 
corrections enhance Born cross sections as shown in Fig.4(a-b), and both the atree and the 
EW corrected atotal in Fig-4(a) are at least one order larger than in Fig.3(a). atotal can 
reach maximal value 0.33 fb at c.m.s energy near 700 GeV and decreases smoothly with 
the increment of ^/s. The corresponding relative correction is depicted in Fig. 4(b) in the 
same parameter space, the correction varies gently from 9% to 21.4%. 

In Fig.5(a-b) we depict the Born and EW corrected cross sections and corresponding 
relative corrections as the functions of mass M^o with y/s = 1000 GeV and the relevant 
MSSM parameters adopting from the SPSla' point except Higgs sector for process e^e~ — >■ 
h^H~^W~. The other Higgs boson masses are determined in the MSSM by using FormCalc 
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Figure 4: The Born and full 0{aew) electroweak corrected cross sections as well as the 
corresponding relative corrections Stotal ^ the functions of the c.m.s energy -s/s with M^o = 
200 GeV and related MSSM parameters being taken from the SPSla' scenario except Higgs 
sector for process e'^e~ —>■ h^H~^W~. (a) cross section versus ^/s, (b) Stotal versus ^/s. 

program with M^o and OS tan /3 = 10.31 being its input parameters. The values of M^o 
and Mjj+ corresponding to different M^o values are also scaled on x-axes in Fig.5(a-b). 
In Fig.5(a), we see there is an abrupt change on each curves due to the resonance effect 
of H+ h'^W+ at the position around M^o ~ ISOGeV (where M^o = 116.8 GeV and 
Mh+ = 197.2 GeV). After reaching its maximal value at position of M^o = 210 GeV, 
the Born electroweak corrected cross sections decrease gently. The dependence of relative 
correction on M^o(or Mf^o, Mh+) is displayed in Fig. 5(b). Stotal curve bends downward 
in the vicinity of M^o ^ 180 GeV because of resonance effect, and then decreases as M^o 
goes up. The relative correction varies from positive to negative value region and the 
absolute value exceeds 10% for both light and heavy Higgs bosons. 

We depict the Born and full one-loop electroweak corrected cross sections as the func- 
tions of c.m.s energy ^/s with the MSSM parameters at the SPSla' point in Fig. 6(a) 
and Fig.6(b) for process e+e~ — > H^H^W~ and e+e~ — > A^H^W~ , separately. In the 
SPSla' scenario, the value of heavier neutral and charged Higgs boson masses(M^o, M^o 
and Mjj+) have the values beyond 420 GeV, accordingly the threshold energy of this pa- 
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Figure 5: The Born and full 0{aew) clcctrowcak corrected cross sections as well as the 
corresponding EW relative corrections Stotal SlS the functions of My^o with -^s = 1000 GeV 
and relevant MSSM parameters taken from SPSla' point except Higgs sector being ob- 
tained by running FormCalc program with input parameters: Mjjo and OS tan/3 = 10.31, 
for process e+e~ — h'^H'^W. (a) cross section versus M^o, (b) relative correction versus 
M40. 

rameter choice can reach 950 GeV. Therefore, in Fig6(a-b) we plot the c.m.s. energy y/s 
in the range between 1000 GeV and 1800 GeV. We see both of Born and 0{aew) corrected 
cross sections increase with the increment of the colliding energy. Fig.6(c) and Fig.6(d) 
show the dependence of the relative corrections 5total of processes e^e~ H'^H'^W' and 
e'^e~ — > A'^H'^W on the c.m.s energy i/s respectively, with the same MSSM parameters 
as taken in Fig.6(a-b). There we see that when the colliding energy -^s goes up, both 
curves in Fig.6(c-d) for relative correction dtotal increase rapidly in the vicinity of thresh- 
old energy, and then becomes increasing smoothly with Stotal value being in the range of 
[—16%,— 23%] when c.m.s energy is beyond 1.2 TeV. 

The Born and full 0{aew) corrected cross sections as the functions of c.m.s energy 
■^s with M^o = 200 GeV and relevant SPSla' MSSM parameters except Higgs sector for 
processes e+e~ H^H+W~ and e+e~ A^H+W~ are plotted in Fig.7(a) and Fig.7(b), 
separately. Again, the values of other Higgs boson masses are obtained from running 
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Figure 6: The Born and full 0{aew) electrowcak corrected cross sections as well as the 
corresponding relative corrections dtotal ^ the functions of the c.m.s energy ^/s in the 
SPSla' scenario for process e+e~ — H'^{A'^)H'^W~ . (a) cross section versus -y/i for 
e+e~ H^H~^W~, (b) cross section versus ^/s for e+e" A^H~^W~, (c) Stotal versus 
for e+e- ^ H^H+W, {d)dtotal versus ^/s for e+e" ^ ^0/^+1^-. 
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FormCalc program with input parameter values of M^o = 200 GeV and OS tan/3 = 10.31. 
Both Born cross section and corrected section increase gently when -y/i goes beyond the 
threshold, while when i/s approaches to 1.2 TeV, the 0{aew) EW corrected cross section 
can achieve the maximal value of about 0.4 fb. The corresponding EW relative corrections 
with the same conditions as in Fig.7(a) and Fig.7(b), are described in Fig.7(c) and (d), 
respectively. As we can see in these figures, the line-shapes of Born and EW corrected 
cross sections (relative correction) for e'^e" — H^H'^^W" process are very similar with the 
corresponding ones for e'^e" A^H''^W~ process. 

The Born and EW corrected cross sections as the functions of input parameter M^o 
are drawn in Fig.8(a) for process e+e~ — *■ H^H~^W~ , and the corresponding relative cor- 
rections versus M^o are depicted in Fig.8(b), with c.m.s energy -s/s = 1200 GeV and 
the relative MSSM parameters taken from SPSla' point except the Higgs sector be- 
ing fixed by using FormCalc program. In Fig.8(a) the one-loop electroweak corrected 
cross section decreases from 0.56 fb to 0.05 fb with the increment of Mjj+{Mf^o) from 
170.28 GeV{152.m GeV) to 408 Gey(400.51 GeV) for process e+e" ^ H^H+W. 
Fig.8(b) shows the relative correction for process e+e~ —>■ H'^H~^W~ goes from —7.8% to 
-20.6% when Mh+{Mho) goes up from 170.28 Gey(152.96 GeV) to 408 GeF (400.51 GeV). 
For channel e'^e" — ^ H~^W~' , the cross section and EW relative correction versus M^o 
are depicted in Fig.8(c,d) respectively, with ^/s being 1000 GeV and other MSSM param- 
eters taken from SPSla' point except Higgs sector being determined by using FormCalc 
package. We can read from Fig.8(c-d) that the EW corrected cross section goes down from 
0.56 fb to 0.005 fb and the relative correction varies from —7% to —28.6%, when input 
value of M^o increases from 150 GeV to 400 GeV. In these four figures of Fig.8, we can 
see the behavior of these two processes e~^e~ H^[A^)H~^W~ are very similar except 
tiny different in the numerical value. 

From above numerical results we can see that the experimental statistic errors for 
processes e+e" h^H+W~, e+e~ H^H+W and e+e~ A'^H+W~ can reach 5.5%, 
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Figure 7: The Born and full O^a^w) clectroweak corrected cross sections as well as the 
corresponding relative corrections dtotal as the functions of the c.m.s energy ^/s with M^o = 
200 GeV and relevant SPSla' MSSM parameters except Higgs sector for processes e+e~ — 
H^{A^)H~^W^ .{r) cross section versus -s/s for e~^e~' H^W^ , (b) cross section versus 

^/s for e+e^ A^H^W~, (c) Stotal versus -s/s for e+e" — >■ H^H~^W~, {d)6total versus -s/s 
for e+e- A^H+W'. 
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Figure 8: The Born and full 0{aew) electroweak corrected cross sections as well as the 

corresponding relative corrections dtotal as the functions of iW^o with c.m.s energy ^/s = 

1200 GeV and relevant MSSM parameters in SPSla' except Higgs sector, for processes 

e+e" H^{A^)H-^W~. (a) cross section versus My^o for e+e" H^H+W~, (b) Stotal 
- uO u+u/- + o- 40 IJ+U/- 



4.2% and 4.2%, respectively. If we assume = 0.6 - 1.8 TeV, J Ldt ~ 1000 
and do not consider the backgrounds to the e+e" = h^^H^^A?) signals 

and the systematic errors of the experiment (due to lacking their studies at ILC), the EW 
corrections to the cross sections of these three processes should be taken into account at 
a linear collider with realistic accuracy better than 10%. 

4 Summary 

In this paper, we present the calculation of the full 0{af>yj) electroweak corrections to 
the processes e+e" ^^H^W^{^ = h^,H^,A^) at a LC in the MSSM. In the nu- 
merical calculation, we generally adopt the MSSM parameters at the reference point 
SPSla' defined in the SPA project. For channel e+e~ h'^H'^W^, we present the 
details for treating the charged Higgs boson resonance in the calculation up to one- 
loop level. Our numerical results show that the EW relative correction for the process 
e+e" — > h^H^W^ can be either positive or negative with the increment of Higgs boson 
mass -M40, and the relative correction varies in the range of [—15%, 20%] as the M^o 
goes from 150 GeV to 400 GeV. For channels e+e~ H'^{A^)H^W'^ , the EW correc- 
tion generally reduce the Born cross sections and the relative correction is typically of 
order —10% ~ —20%. We conclude that the complete 0{aew) electroweak corrections to 
the process e+e~ —>■ ^^H'^W^{^ = h^,H^,A^) are generally significant and cannot be 
neglected in the precise experiment analysis. 
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